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I.  INTRODUCTION 


Disturbances  are  defined  as  the  uncontrollable  inputs  which  act  on  a  dyna¬ 
mical  system.  There  are  many  varieties  of  disturbance  inputs  which  can  be 
associated  with  a  controlled  system  and  they  are,  for  the  most  part,  comple¬ 
tely  unpredictable  in  magnitude  and  in  their  arrival  times. 

In  practice,  additive  disturbances,  i.e.,  disturbances  which  are  repre¬ 
sented  by  terms  added  to  the  plant  state  equation,  can  arise  from  motivating 
effects  external  to  the  plant  (external  disturbances)  or  from  motivating 
effects  arising  from  the  physical  characteristics  of  plant  subsystems  or 
internal  plant  dynamics  (internal  disturbances).  Further,  these  disturbances 
can  be  divided  into  two  categories:  (a)  noise  disturbances,  characterized  by 
random  and  erratic  behavior  with  relatively  high-frequency  content,  and  (b) 
waveform  structured  disturbances,  characterized  by  a  degree  of  waveform  regu¬ 
larity  which  can  be  described,  piecewise  in  time,  by  differential  equations 
forced  by  sparse  sequences  of  impulses.  The  nature  of  these  disturbances  may 
be  either  completely  known  (through  direct  prior  or  real-time  observation  or 
test),  completely  unknown  (random-like),  or  partially  known. 

Johnson  [1-6,  10]  introduced  the  idea  of  mathematically  describing  uncer¬ 
tain  waveform-structured  disturbances  by  representing  them  as  a  weighted 
linear  combination  of  known  basis  functions  of  the  form 

n 

w(t)  -  S  cJLf1(t),  (1) 

i-1 

where  w(t)  is  the  plant  disturbance  vector  and  is  a  p-vector  and  the  weighting 
coefficients  c^  are  completely  unknown  constants  which  can  change  in  magnitude 
in  a  random,  once-in-a-while  fashion.  The  basis  functions  f^(t)  are  complete¬ 
ly  known  because  they  are  chosen  by  the  designer  based  on  the  waveform  pat¬ 
terns  exhibited  (or  thought  to  be  exhibited)  by  the  disturbance. 

Johnson  [1-11]  developed  a  control  engineering  design  technique,  referred 
to  as  disturbance  accommodation,  wherein  a  combination  of  waveform-mode 
disturbance  modeling  and  state-variable  control  methods  are  utilized  to  design 
controllers  which  will:  (1)  absorb  (counteract),  (2)  minimize,  or  (3)  construc¬ 
tively  utilize  the  effects  of  uncertain  disturbances  on  the  plant.  Three  main 
classes  of  controllers  are  considered  within  the  overall  cognomen  of  distur¬ 
bance  accommodating  control  theory.  These  are:  (1)  Disturbance  Absorption 
Controllers  (DAC),  (2)  Disturbance  Minimization  Controllers  (DMC),  and  (3) 
Disturbance  Utilizing  Controllers  (DUC).  Each  class  of  controller  has  its  own 
associated  design  goals  and  design  methodology.  The  mathematical  theories  of 
DAC  and  DUC  were  thoroughly  developed  in  References  1  through  12.  The  theory 
and  techniques  associated  with  DMC  were  compiled  and  extended  in  Reference  13. 
Additional  results  pertaining  to  the  application  of  DMC  techniques  to  a 
linear,  time-invariant,  second-order  state  set-point  regulator  problem  with 
constant  external  disturbances  were  presented  in  Reference  14. 

One  of  the  DMC  techniques  presented  in  Reference  13  was  the  isobasls 
design  technique.  This  technique  makes  the  a  priori  assumption  that  the 
control  vectors  will  be  some  function  of  the  disturbance  vectors,  i.e.,  they 
will  be  composed  of  some  combination  of  the  same  basis  functions  which 
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describe  Che  disturbances.  This  approach  provides  a  parametric  form  for  the 
control  vectors.  By  utilizing  the  preferred  minimization  method  [13],  one  can 
design  the  parameters  of  the  minimization  controller  to  minimize  the  distur¬ 
bance  effects. 

Two  examples  of  the  application  of  the  Isobasis  technique  were  presented 
in  Reference  13.  One  example  was  for  a  second-order  output  servo-command 
prob.lem  with  a  ramp  for  the  external  disturbance  and  a  ramp  for  the  output 
servo-command.  The  other  example  was  for  a  second-order  state  servo-command 
problem,  again  with  a  ramp  for  the  external  disturbance  and  for  the  inpuc 
state  servo-command.  This  report  will  present  results  on  the  application  of 
the  Isobasls  design  technique  to  linear,  time-invariant  state  set-point  regu¬ 
lators  with  time-varying  external  disturbances. 

II.  LINEAR  DYNAMICAL  SYSTEMS 

The  class  of  systems  to  be  considered  in  this  report  are  "linear, 
time-invariant,  dynamical  systems,”  so  called  because  the  vector  differential 
equation  for  the  state  x(t)  is  a  linear  differential  equation,  the  transfor¬ 
mation  between  the  state  space  and  output  space  is  linear,  and  the  elements  of 
the  matrices  in  the  plant  model  are  constant  with  respect  to  time. 

These  systems  will  be  represented  by  equations  of  the  general  form 

x(t)  ■  Ax(t)  +  Bu(t)  +  Fw( t)  (2) 

y(t)  -  Cx(t)  +  Eu( t)  +  Gw(t)  (3) 

where  x(t)  is  the  plant  state  vector  and  is  an  n-vector,  u(t)  is  the  plant 
control  input  vector  and  is  an  r-vector,  w(t)  is  the  plant  disturbance  vector 
and  is  a  p-vector,  y(t)  is  the  plant  output  vector  and  is  an  m-vector  and  A, 

B,  F,  C,  E,  and  G  are  appropriate  size,  known  matrices  with  time-invariant  el¬ 
ements.  In  addition,  the  general  form  of  the  disturbance  state  model  is  [10]. 

w(t)  -  Hz(t)  +  Lx(t)  (4) 

z(t)  ■  Dz(t)  +  Mx(t)  +  a(t)  (5) 

where  z(t)  is  the  p-dimensional  disturbance  state  vector,  cr(t)  is  a  sparsely 
populated  vector  impulse  sequence  and  H,  L,  D,  and  M  are  appropriate  size, 
known  matrices. 
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III.  BACKGROUND 


In  Reference  13,  several  methods  were  presented  for  minimizing,  via 
direct  control  action,  the  effects  of  constant  disturbance  components  which 
are  not  completely  absorbable  on  linear,  time-invariant  state  set-point 
regulators.  The  metric  used  for  the  minimization  process  is  the  norm  defined 
by 

|  |  Ax  -  b  |  |p  -  (Ax-b)T  Q(Ax-b) ,  Q  >  0.  (6) 

The  design  objective  is  the  minimization  of  the  distance  between  the 
attainable  and  desired  set-point,  where  this  distance  is  defined  by  the 
Euclidean  norm, 

d2  ,||  £  ||2  .  eTe  f  (7) 

of  the  error  vector  between  xari  and  the  plant  state  x(t),  i.e., 

e(t)  -  x8p  -  x(t).  (8) 

An  expression  for  the  dynamics  associated  with  this  error  can  be  derived  by 
differentiating  Equation  (8)  and  substituting  in  the  appropriate  terms  from 
Equation  (2).  The  result  can  be  expressed  as 

*e(t)  -  xgp  -  x(t)  -  Ae(t)  -  Bu(t)  -  Axsp  -  Fw(t)  ,  (9) 

where  Axgp  represents  the  "set-point  disturbance"  term. 

In  disturbance  accommodating  control  design,  the  control  vector  u(t)  is 
considered  to  be  an  ordered  collection  of  the  various  independent  control 
Inputs  which  are  available  to  accomplish  the  primary  control  objective  and  to 
“accommodate"  the  disturbances  which  are  acting  on  the  system.  In  the  design 
of  disturbance  minimization  controllers,  it  is  common  practice  to  split 
(allocate)  the  total  control  u(t)  into  two  parts  as  follows, 

u(t)  -  up(t)  +  ud(t)  ,  (10) 

where  up(t)  is  given  the  task  of  accomplishing  the  primary  control  objective 
and  ud(t)  is  given  the  task  of  disturbance  accommodation.  The  part  ud(t)  can 
be  further  allocated  into  component  vectors,  as  required.  For  the  methods 
considered  in  this  report,  ud(t)  will  sometimes  be  allocated  as 

ud(t)  -  ud8(t)  +  udw(t)  *  UD 

The  component  uds(t)  will  be  designed  to  accommodate  the  effects  of  the  set- 
point  disturbance  term,  while  udw(t)  will  be  designed  to  accommodate  the 
effects  of  the  external  disturbance  term.  If  the  plant  is  completely 
controllable  and  is  also  completely  observable,  the  control  up(t)  can  be 
designed  in  the  form 


Given  the  allocation  of  the  control  vector  u(t),  Equation  9  can  be 
rewritten  as 

e(t)  -  Ae(t)  -  BUp(t)  -  Buds(t)  -  Budw(t)  -  Axsp  -  Fw(t).  (13) 

Upon  substitution  from  Equation  (4),  Equation  (13)  becomes 
*e(t)  -  Ae(t)  -  Bu_(t)  -  Buds(t)  -  Budw(t) 

(14) 

-  Axgp  -  FHz(t)  -  FLx8p  -  FL e( t) . 

In  the  case  of  Equation  (14),  one  would  design  up(t)  in  the  form 

up(t)  -  -  Ke(t)  ,  (15) 

with  K  chosen  such  that  the  homogeneous  system 

*e(t)  -  (A  +  FL  +  BK)e(t)  (16) 

will  yield  e(t)  — 9  0  "rapidly."  If  one  lets  A  ■  A+FL+BK,  then  Equation  (14) 
can  be  expressed  as 


e(t)  -  Ae(t)  -  [(A-fFL)xSp  +  Buds(t)J  -  [FHz(t)  +  Budw(t)] 


(17) 


IV.  PLANT  AND  EXTERNAL  DISTURBANCE  MODELS 


The  plant  state  and  output  models  to  be  used  in  the  examples  presented  in 
this  report  are 


CH  *CK)' 


The  target  state  set-point  vector  is  given  as 


(18) 

(19) 


xsp  *  ^*sp,l*  (20) 

The  plant  given  by  Equations  (18)  and  (19)  is  completely  controllable.  Por 
the  purposes  of  the  examples  it  is  assumed  that  all  necessary  state  infor¬ 
mation  is  available  from  an  ideal  state  reconstructor. 


Several  different  types  of  external  disturbance  will  be  applied  to  the 
plant  and  the  effectiveness  of  the  disturbance  minimization  controllers  on 
each  will  be  examined.  The  disturbances  which  will  be  used  and  the 
corresponding  models,  based  on  Equations  (4)  and  (5)  are  as  follows. 


1.  w(t) 


w 


3.  w(t) 


w 


C0  +  Ctt 
(1,  0)z 


(1,  0)z 
0  1 
-a?  0 


+  cr(t) 


+  c(t) 


C0  +  cie<* 

(1,  0)z 


(21) 

(22) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


V.  GEOMETRICAL  CONSIDERATIONS  AND  THEOREMS 


Since  Che  control  distribution  matrix  B  of  Equation  (18)  is  a  2x1  matrix 
of  rank  1,  it  does  not  span  the  state  space,  which  is  two-dimensional  in  this 
example.  Hence,  Bud(t)  will  have  a  limited  set  of  attainable  points  in  the 
state  space.  Also,  the  external  disturbance  distribution  matrix  is  of  rank  1 
and  thus,  Fw(t)  will  have  a  limited  range  of  action  in  the  state  space.  As 
can  be  seen  in  Figure  1  the  lines  of  action  of  Bud  and  Fw  are  not  colinear. 
Hence,  no  ud  exists  which  will  completely  absorb  a  non-zero  external  distur¬ 
bance  or  a  set-point  disturbance  resulting  from  a  non-zero  target  state  set- 
point. 

Given  that  this  situation  exists  and  that  a  design  objective  is  to  mini¬ 
mize  the  effects  of  the  uncancellable  disturbances,  one  thus  attempts  to 
design  ud  so  as  to  achieve  this  objective  in  some  fashion.  With  respect  to 
the  vectors  Fw^  and  Budi  shown  on  Figure  1,  one  approach  to  the  minimization 
problem  is  to  first  express  the  vector  Fw^  as  the  sum  of  two  component  vectors, 
one  lying  in  the  column  range  space  of  B,  R(B),  and  one  lying  in  the  orthogo¬ 
nal  complement  to  the  column  range  space  of  B,  R(B)X.  This  makes  it  easy  to 
see  that  the  component  lying  in  R(B)1,  which  is  the  component  that  is 
uncancellable,  is  minimized  if  the  component  lying  in  R(B)  is  the  orthogonal 
projection  of  Fw^  onto  R(B)  [13 J.  How  can  udi  be  chosen  such  that  this  is 
accomplished? 

Casting  the  problem  into  the  form  of  Equation  (6)  one  wishes  to  minimize 

| | Budl  +  Fwl  ||q  “  <Budl  +  Fw1>T  Q<Budl  +  Fwl>  •  (3°) 


Consider  the  following  theorems  from  [15]. 


Theorem:  Let  B  be  an  mxn  matrix,  Fw^  an  m  vector  an  Q  a  positive  definite 
mxm  matrix.  Then  |  j Budi  +  Fw^  |  |  is  smallest  when 

udl  -  XFwx  ,  (31) 

where  X  satisfies 

BXB  -  B  ,  (QBX)T  -  QBX  .  j  (32) 

Theorem:  Let  B  be  an  mxn  matrix,  Fw^  an  m  vector  and  P  a  positive  defi¬ 
nite  nxn  matrix.  If  Budi  ■  -Fw^  has  a  solution  for  udi»  the  unique  solution 
for  which  ||udi  | j  p  Is  smallest  is  given  by 

udi  -  XFWl  ,  (33) 

where  X  satisfies 


BXB  -  B  ,  (PXB)T  -  PXB . 


(34) 
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Corollary:  Let  B  be  an  mxn  matrix,  Fwj  an  m  vector,  Q  a  positive  definite 
axm  matrix  and  F  a  positive  definite  nxn  matrix.  Then,  there  is  a  unique 
matrix  X  satisfying 

(QBX)T  -  QBX  ,  (PXB)T  -  PXB  .  (35) 

Moreover,  | | Budi  +  Fw^  II  q  assumes  its  minimum  value  for  ud^  *  XFw^,  and  in 
the  set  of  vectors  ud^  For  which  the  minimum  value  is  assumed,  udi  “  XFw^  is 
the  one  for  which  |  |  ud^  | | p  is  smallest.  | 

If  the  weighting  matrices  Q  and  P  are  chosen  to  be  the  identity  matrix  I, 
then  Equation  (30)  represents  the  Euclidean  norm  and  the  properties  which  X 
satisfies  in  the  above  theorems  are  properties  which  are  satisfied  by  the 
Moore-Penrose  generalized  inverse.  Thus,  if  udj^  is  chosen  as 

udl  -  -B^Fwj  ,  (36) 

where  ( •  V  represents  the  Moore-Penrose  generalized  inverse  of  (O,  then 

udi  is  the  minimum  norm  controller  which  minimizes  |  j Bud^  +  Fw^  ||  .  Also,  one 

has 

Budl  +  Fw1  -  -BB^FW!  +  Fwx  -  (I  -  BBt)Fw1  (37) 

and  (I-BB^)  is  the  projector  of  Fw^  on  R(B)1  along  R(B)  and  BB^  is  the  projec¬ 
tor  of  Fwi  on  R(B)  along  R(B)i>.  The  uncancellable  part  of  Fw^  is  thus  the 
component  in  R(B)X. 


VI.  DISTURBANCE  MINIMIZATION  WITH  w(t)  -  c0  +  cxt 

A.  With  Allocated  Disturbance  Minimizing  Control  Vector 

This  section  will  apply  the  isobasis  design  technique  to  a  state  set¬ 
point  regulator  problem  where  the  external  disturbance  is  given  by  Equation 
(21),  l.e.,  a  combination  of  a  step  and  a  ramp.  The  stabilization  control 
Up,  designed  according  to  Equation  (15)  and  Equation  (16),  is  chosen  to  be 


113,  14] 

up(t)  -  -Ke(t)  -  (3.,  0.36) e(t) .  (38) 

Recall  Equation  (17),  which  is  the  general  expression  for  e(t)» 

*e(t)  -  Ae(t)  -  (Ax8p  +  Buds)  -  [Fw(t)  +  Budw(t)].  (39) 

If  w(t)  »  0,  a  steady-state  solution  would  exist  for  e  as  [14] 

%s  *  A-ifAXgp  +  Buda] .  (40) 

This  steady-state  error  could  then  be  minimized  by  choosing  uds  as 

uds  "  <A“1B)tA-lAx8p  (41) 

which,  in  this  example,  becomes  [13,14] 

uds  "  -2*3447xap,l  •  («) 


In  order  to  develop  an  expression  for  u^t),  with  the  external  distur¬ 
bance  given  by  Equation  (21),  the  a  priori  assumption  is  made  that  udw  will  be 
of  the  form 


udw(t)  “  b<co  +  clfc)  "  ^(t)* 

The  general  solution  of  Equation  (17)  for  e(t)  can  be  written  as 
e(t)  -  eAte(0)  -  f  eA(t-t)[Ax8p  +  Buds]dx 

*4> 

/t  ~ 

eA(t-x)  [Fw(  x)  +  Budw(x)]dx. 


(43) 


(44) 


If  Equation  (21)  and  Equation  (43)  are  substituted  into  the  last  term  on  the 
righthand  side  of  Equation  (44),  the  error  contribution  due  Co  the  external 
disturbance  term  can  be  expressed  as 


/t  ~  r  t  ~ 

eA(t-x)  (f  +  Bb)c0dT  -  J  eMt-x)  (p  +  Bb)ctxdT 
o  o 

A#  M 

-  -  A-l(eAt-i)(F+Bb)c0  -  [A“2(eAt-i)  -  A“lt] (F+Bb)cx  .  (45) 

If  one  considers  the  response  after  the  transients  have  settled  out.  Equation 
(45)  can  be  reexpressed  as 

A^  -  A-1(F  +  Bb)c0  +  A“1(A“1  +  It)(F  +  Bb)c1  .  (46) 

Upon  expanding  the  right-hand  side  of  Equation  (46),  one  obtains 


A£W 


0854  -0.1951 

829  -0.6098 


[:: 

E-0.35  +  0.< 

-0.9591  +  l.J 


(47) 


- 


.0854t 
.829t 

which  can  be  rearranged  to  the  form 
^-0.1097w  -  0.2477w 
1.2l92w  -  0.9441w 
Let  Equation  (48)  be  represented  as 
Ae^  •  u^b  +  U2  • 

One  would  like  to  have 
u^b  +  U2  ■  0 


0.1023  -  0.1951t 
0.015  -  0.6098t 


^-0.3048w  -  0.l45w 
0.6094w  -  0.9291w> 


(48) 


(49) 


(50) 


but  this  is  not  attainable  in  the  case  of  this  example.  If  the  norm  minimiza¬ 
tion  criterion  is  applied,  one  chooses  b  such  that  nAew||is  minimized. 

Doing  this,  that  b  of  minimum  norm  which  will  minimize  the  norm  of  A^  (see 
Section  V)  is  found  to  be 


b*  — (ui)TU2 
where [13] 


u. 


(u,^) 


-1  T 
*•11 


(51) 


(52) 


When  the  appropriate  substitutions  are  made  from  Equation  (48)  into  Equation 
(52)  and  Equation  (51),  the  results  obtained  are  as  follows, 


10 


J  J  *  f 


£ 


.  >  W  "  ^  M  L.  .  ^  ^  L^-  ^  -  •>  - 


r« 


t 

U1 


b* 


[(-0.3048w  -  0.1454w,  0.6098w  -  0.9291w)  ("oieo^  -  o!929l£)]'1Ul 

-0.3Q48w  -  0.l454w  ,  _ 0.6094w  -  0.9291w 

4643w2  -  1.04 38ww  +  0.8843w^  0.4643W2  -  1.0438ww  +  0.8843w2 

-0.7764w2  +  1.6167ww  -  0.9132w2 
0.4643w2  -  1.0438ww  +  0.8843W2 


(= 


)• 


(53) 

(54) 


For  the  purposes  of  these  examples,  it  will  be  assumed  that  w  and  w  are 
obtained  from  an  ideal  disturbance  state  reconstructor. 

A  digital  simulation  was  written  for  this  example.  Several  runs  were  made 
with  various  values  for  the  coefficients  describing  w(t).  The  results  are 
shown  in  Figures  2  through  9.  In  each  case,  three  curves  are  given.  One 
curve  represents  a  case  where  u^  ■  0.  One  curve  is  for  the  controller 
designed  in  this  section.  A  third  curve,  which  is  shown  for  purposes  of 
comparison,  is  for  a  case  where  u^w  was  designed  under  the  assumption  that  w 
was  a  constant  disturbance.  With  this  assumption,  when  u<iw  is  designed  to 
minimize  the  contribution  of  the  external  disturbance  to  the  steady-state 
error,  the  result  is 


u^n  *  — 1.6723z. 


(55) 


Figures  2  through  5  show  the  resulting  set-point  error  magnitude  for 
various  target  set-points  and  external  disturbances.  Figures  6  through  8  show 
the  set-point  error  magnitude  and  the  error  components  for  a  case  with  zero 
target  set-point  and  changing  external  disturbance  parameters  and  Figure  9 
shows  w(t).  As  can  be  seen  from  the  plots,  the  controller  with  udw  designed 
using  the  lsobasls  technique  with  norm  minimization  does  not  result  in 
improved  performance  over  the  controller  designed  to  minimize  e88  when  w  is 
assumed  constant. 

Another  minimization  criterion  which  could  be  applied  would  be  to  design 
u<iw  such  that  A%,i  or  Aew2  is  steered  to  zero,  l.e.,  design  u^w  such  that  the 
direct  effect  of  the  external  disturbance  on  ei  or  62*  respectively,  is 
completely  absorbed.  If  one  follows  this  approach,  the  parameter  b  would  be 
designed  as  follows.  First,  assume  that  ei  is  a  critical-state  variable  and 
that  the  direct  effect  of  w  on  ei  is  to  be  absorbed.  From  Equation  (48),  one 
has  that 


A«wl  -  -(0.1097w  +  0.2477w)  -  (0.3048w  +  0.1454w)b. 


In  order  to  obtain  Aewi  -  0,  one  thus  requires  that 
b  1 


0.1097w  +  0.2477w 
0.3048w  +  0.1454w 


(56) 


(57) 
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Figure  3.  Error  versus  time,  with  allocated  control  vector, 
w-l+t,  Xsn-  (10.0)1. 
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Figure  8.  £2(t)»  case  with  parameter  variation  in  w(t) 
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Figure  9.  External  disturbance  time  history. 


«.  V  V->.  ' 


l 


i 


I 

i  * 

l>  ' 

! 

S 


Figures  10  through  15  show  results  for  two  cases  when  Equation  (57)  is  used  to 
define  udw.  As  can  be  seen  in  Figures  11  and  14,  in  each  case  e^(t)  was 
steered  to  zero.  However,  the  total  set-point  error  (Figures  10  and  13)  was 
not  Improved  over  that  obtained  by  using  u^  designed  to  minimize  e88. 

Next,  assume  that  *s  the  critical-state  variable  and  design  u<jw  to 
absorb  the  direct  effect  of  w  on  62 •  Again,  from  Equation  (48),  one  has 

A €^2  -  (1.2191w  -  0.944lw)  +  (0.6094w  -  0.929lw)b.  (58) 

In  order  to  obtain  Ae^j  "  0»  one  requires  that 
*  _  1.2192w  -  0.944lw 

d  *  -  ———————  .  (jy) 

0.6094w  -  0.9291w 


Figure  16  shows  results  obtained  for  the  case  used  in  Figures  10  through  13 
when  Equation  (59)  is  used  to  define  u<jw.  As  can  be  seen,  this  controller 
produced  unacceptable  transients.  This  points  up  a  caution  in  the  use  of  this 
technique. 

Since  the  parameters  given  by  Equation  (57)  and  Equation  (59)  to  define 
udw(t)  have  time  varying  functions  in  their  denominators,  the  possibility 
exists  that  some  comblnation(s)  of  the  constants  cQ,  c^  defining  w(t)  may 
result  in  a  zero  in  the  denominator  thus  causing  an  infinitely  large  value  to 
be  input  by  the  disturbance  minimizing  controller.  For  the  two  external 
disturbances  used  to  generate  the  plots  of  Figures  10  through  16,  i.e.. 


w(t)  ■  1  +  t 

(60) 

w(t)  *  -2  -5t 

(61) 

a  check  of  the  denominator  of  Equation  (57)  shows  that  no  positive  value  of 
time  exists  at  which  u<jw  would  go  to  infinity.  However,  the  denominator  of 
Equation  (59)  will  go  to  zero  at  t  ■  0.525  seconds  for  Equation  (60)  (which 
gives  the  resulting  peak  in  the  response  shown  in  Figure  16)  and  at  t  -  1.125 
seconds  for  Equation  (61).  A  practical  implementation  would  thus  require  a 
limit  on  the  allowable  magnitude  of  u,jw.  Figures  17  through  19  show  the 
results  for  w  as  in  Equation  (60)  when  the  magnitude  of  u<jw  is  limited  to  +50. 
From  Figure  19  it  can  be  seen  that  62^)  is  steered  to  near  zero  and  a  com¬ 
parison  of  Figures  17  and  10  indicates  that  choosing  u<iw  to  give  Ae^  ■  0  does 
reduce  the  total  error  over  the  case  where  Ujjw  is  chosen  to  give  Ae^i  “  0. 
However,  in  neither  case  is  the  total  error  less  than  that  obtained  using  the 
controller  designed  to  minimize  Cgg. 

B.  With  Unallocated  Disturbance  Minimizing  Control  Vector 

In  Fart  A,  the  disturbance  minimizing  control  u<j  was  allocated  as  u<j  “ 
uds  +  udw  aQd  ud8  wa8  designed  to  minimize  the  norm  of  the  steady-state  set- 
point  error  contribution  given  by  Equation  (40).  As  was  seen,  the  performance 
of  the  controllers  designed  in  Part  A  was  not  as  good  as  that  of  the  control¬ 
ler  designed  to  minimize  the  steady-state  error  under  the  assumption  that  w(t) 
was  a  constant.  In  this  section,  the  isobasis  design  technique  will  be 
applied  with  an  unallocated  disturbance  control  vector. 
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Figure  10.  Error  versus  time,  udw  designed  to  give  Aewi=0. 


2 . 00 


4.00 


.  ©0 


8.00 


10.00  12.00 


T IMECSECD 

11.  £l  (t)  with  uju  designed  to  give  Ae„,-0. 


3  2.00 

T  IMEC 


4.00  6.00  8.00  J0.O0  12.00 

■  EC  3 


av,  designed  to  ftive  -  0. 


Figure  13.  Error  magnitude,  u, 
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Figure  16.  Error  versus  time,  u-dw  designed  to  give  AeM?=0., 
no  limit  on  ujv. 


Id  4.(00  6.00  8.00  10.00  12.00 

TIMECSECD 


Figure  17.  Error  versus  time,  udw  designed  to  give  Aew2=0., 
udw  limited  to  ±50. 


Since  the  set-point  error  is  a  constant,  the  a  priori  assumption  will  be 
made  that  ud  is  of  the  form 

ud(t)  -  b(cQ  +  ^t  +  c2),  (62) 

where  c2  is  included  to  represent  the  fact  that  the  set-point  disturbance  is 
constant.  The  general  expression  for  e(t)  is  thus  given  as 

e(t)  -  Ae(t)  -  [Ax3p  +  Fw(t)]  -  Bud(t).  (63) 

The  general  solution  of  Equation  (63)  for  e(t)  can  be  written  as 


e(t) 


eAte(O) 


IAxsp  +  Fw(  r)  +  Bud('t)]dx 


eAte(0)  -  J  £A(t-T)[Axsp  +  Fcq  +  Bb(cQ  +  c2)]dt 


-/ 


t  ~ 

eA(t-x)  (Fci  +  Bbc^Tdt 


(64) 


eAtE(0)  -  A“l(eAt  -  I)[Axap  +  FcQ  +  Bb(cQ  +  c2)] 
-  [A-2(eAt  -  I)  -  A_1t] (F  +  Bb)c1 


Again,  if  one  considers  the  response  after  the  transients  have  settled 
out.  Equation  (64)  can  be  reexpressed  as 

e(t)  -  A-ifAXgp  +  FcQ  +  Bb(c0  +  c2)]  +  A-UA-1  +  It)(F  +  Bb)ci  .  (65) 


When  the  appropriate  substitutions  are  made  into  Equation  (65),  the  result  is 
0.0854  -0.1951" 


e(t) 


1.829 


-0.6098 


rxSp,l  +  bc2  +  (l+b)c0\ 
2bc2  +  (l+2b)c0  / 


-0.35  +  0.0854t  0.1023  -  0.1951t 

-0.9591  +  1.829t  0.015  -  0.6098t  . 


1+b 


a+2by 


Cl 


'0.0854xsp>l  -  0.1097w  -  0.2477w 
^1.829xgp>i  +  1.2l92w  -  0.9441w 


(66) 


+ 


-0.3048C2  -  0.3048w  -  0.1454w 
.  0.6094c2  +  0.6094w  -  0.9291w 


b. 


If  Equation  (66)  is  written  in  the  form 


e(t)  •  u^b  +  U2  (67) 

then  that  b  of  minimum  norm  which  will  minimize  the  norm  of  e(t)  is  found  to 
be  given  by 

b*  ■  -(u^)^U2  “  -(u^ui)-^  u^u2  •  (68) 


When  Equation  (68)  is  evaluated,  the  resulting  expression  for  b*  is, 

b*  - 

[1.0886(c2-»w)-1.7117w]x3Ptl.K)t7764w2.).(0t7764w_0t4498;;)e2,lt6166w^K)>9132y2] 

0 . 4643c2+0 • 9286C2W-1 . 0438c2w+0 • A643w2-1 . 0438ww+0 . 8843W2 

(69) 

It  should  be  noted  that  if  w(t)  is  assumed  to  be  zero,  the  corresponding 
controller  using  b*  from  Equation  (69)  is 

ud  -  -2.3447xap>i  ,  (70) 

which  is  the  result  obtained  in  Equation  (42). 

The  digital  simulation  from  Part  A  was  used  to  generate  results  with  b*  as 
given  by  Equation  (69).  The  two  cases  given  by  ud  *  0  and  ud  designed  to 
minimize  EgS  were  again  Included  for  purposes  of  comparison.  Results  are  pre¬ 
sented  in  Figures  20  through  24. 

As  shown  in  Figure  20,  when  w(t)  is  a  constant  the  controller  from  this 
section  and  that  designed  to  minimize  ess  give  identical  results.  Figures  21 
through  24  repeat  the  cases  shown  in  Figures  2  through  5,  respectively.  A 
comparison  of  Figures  21  and  2  shows  that  with  no  interaction  possible  between 
the  set-point  and  external  disturbances  (since  the  set-point  is  the  origin), 
the  disturbance  minimizing  controllers  of  Part  A  and  Part  B  give  identical 
results.  When  interactions  are  possible,  l.e.,  neither  the  set-point  nor 
external  disturbances  are  zero,  a  diversity  of  results  occurs.  Comparing 
Figures  22  and  3,  it  can  be  seen  that  the  controller  of  Part  B  gave  better 
overall  results  than  did  the  controller  of  Part  A.  A  comparison  of  Figures  23 
and  4  shows  that  the  controller  of  Part  B  resulted  in  much  better  performance 
than  did  the  controller  of  Part  A,  except  around  a  time  of  3.5  seconds,  with 
zero  error  at  two  different  times.  However,  in  comparing  Figures  24  and  5  it 
can  be  seen  that  the  controller  of  Part  B  gave  slightly  better  results  between 
a  time  of  2  and  5  seconds,  and  much  worse  results  between  a  time  of  5  and  9 
seconds,  as  compared  to  the  controller  of  Part  A.  It  is  also  apparent  that, 
except  possibly  for  the  case  shown  in  Figure  23,  the  controller  designed  to 
minimize  EgS  still  results  in  the  best  overall  performance. 
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VII.  DISTURBANCE  MINIMIZATION  WITH  w(t)  -  C!  sin  act 


This  secclon  will  apply  the  isobasis  design  technique  to  the  state  set- 
point  regulator  example  with  an  external  disturbance  of  the  form 

w(t)  ■  c^slnat  .  (71) 

For  this  case,  the  a  priori  assumption  is  made  that  u<jw  is  of  the  form 

Udtf(t)  ■  bc^slnoct  *  bw(t),  (72) 

and  u^8  is  as  shown  in  Equation  (42). 

That  part  of  the  solution  of  Equation  (17)  for  e(t)  which  is  due  to  the 
last  term  on  the  right-hand  side  of  Equation  (17)  will  again  be  found,  with 
ur?v  as  given  by  Equation  (72).  In  this  case  then,  one  has 


e^(t-r)  [yw(  T)  +  Bbw(t)]dx 


eA(t“T)(F  +  Bb)c^sinaxd'c 


(73) 


■  -  (A “2  +  <£i)-l  (-Asinat  -  alcosat  +  eAta)(F  +  Bb)c1  . 


The  response  after  the  initial  system  transients  have  settled  out  will  again 
be  considered  in  designing  udw,  thus,  Equation  (73)  is  rewritten  as 

A^,  -  -&-l(-Aw  -  wI)(F  +  Bb),  (74) 

where  ^  ■  (A2  +  o?I).  From  Equation  (74)  one  can  find  that  that  b  which  is 
itself  of  minimum  norm  and  which  minimizes  the  norm  of  Ae^  is  given  by 

b*  -  -$"1(-Aw  -  wDel^-^-Aw  -  wI]F.  (75) 

When  the  appropriate  substitutions  are  made  into  Equation  (75),  the  following 
is  obtained, 


a_1  ,  _i_  r-3-762  +  02  i-i°i  1 

DET  L  -10.72  0.16  +  a2] 


where  DET  *  oA  -  3.602a£  +10.76  , 


«-l(-Aw  -  wI)F  -  ^ 


(1.182  +  1.36a2)w  +  (2.661  -  a2)w\ 

,  o  •*  (78) 

(-13.12  +  5.72a2)w  +  (10.16  -  a2)w/ 


_  .  /  (3.28  +  0.72a2)w  +  (1.56  -  a2)w\ 

A“l(-Aw  -  wI)B  -  —  j  .  )  ,  (79) 

DET  \(-6.56  +  5.44 az)w  +  (10  -  2a*)w  / 


$“l(-Aw-wI)B]'r  -  tC^“1(-Aw-wI)B]T[^-1(-Aw-wI)B3-l[A-l(-Aw-wi)B]  . 


When  Equation  (75)  is  evaluated,  one  has 
cnw2  +  ci2«w  +  cj3w2 
c2i  *2  +  C22W  +  C23W2 

where 


on  -  89.95  -  103.59a2  +  32.1a4 

ci2  -  -187.28  +  145.55a2  -  18.96a4  (82) 

d3  -  105.75  -  34.54a2  +  3a4 

C21  -  53.79  -  66.65 a2  +  30.11a4 

C22  “  -120.97  +  130.73a2  -  23.2a4  (83) 

C23  -  102.43  -  43.12a2  +  5a4 

The  controller  defined  by  Equation  (72)  using  Equation  (81)  to  Equation 
(83)  was  added  to  the  digital  simulation  of  the  set-point  regulator  example  of 
Section  VI.  Figures  25  through  31  show  results  obtained  with  the  114  of  this 
section  as  compared  to  results  with  114  »  0  and  U4  designed  to  minimize  ess 
under  the  assumption  that  w  is  a  constant.  As  can  be  seen  from  these  results, 
the  external  disturbance  minimizing  control  vector  U4  designed  via  the  isoba¬ 
sis  design  technique  does  provide  better  performance,  in  this  case,  than  U4  as 
designed  to  minimize  es8,  except  at  low  frequencies.  An  examination  of  these 
figures  shows  that  for  the  low  frequency  cases  (Figures  25,  26,  30  and  31 ) 
the  controller  designed  to  minimize  eg8  provided  no,  or  at  least  a  small, 
decrease  in  the  peak-to-peak  amplitude  of  the  sinusoidal  external  disturbance 
effect  over  the  case  with  U4  ■  0.  However,  as  the  frequency  of  the  external 
disturbance  term  Increased,  this  controller  resulted  in  an  Increase  in  this 
peak-to-peak  amplitude.  The  controller  designed  in  this  section,  on  the  other 
hand,  generally  resulted  in  a  decrease  in  this  amplitude.  The  exceptions  were 
for  the  two  cases  where  a  ■  1  rad /sec.  As  the  frequency  increased,  this 
decrease  became  substantial. 

In  order  to  obtain  a  better  indication  of  the  performance  obtained  via  use 
of  the  isobasls  technique  for  a  sinusoidal  external  disturbance,  several  sets 
of  runs  were  made.  The  first  set  was  for  a  case  with  x8p>i  *  10.  Figure  32 
shows  the  percent  by  which  the  peak-to-peak  magnitude  of  the  sinusoidal  com¬ 
ponent  of  the  error  was  reduced,  as  a  function  of  the  frequency  of  w(t),  by 
using  U4  as  designed  in  this  section.  Figures  33  and  34  show  the  percent 
reduction,  as  a  function  of  the  amplitude  of  w(t),  for  frequencies  of  3  and  10 
rad/sec,  respectively.  Figures  35  through  37  present  similar  data  for  a  case 
with  Xsp,l  ■  "5.  All  data  in  each  case  was  obtained  after  the  initial  tran¬ 
sients  had  settled  out.  As  can  be  seen  from  examination  of  Figures  32  and  35, 
for  frequencies  above  2  rad/sec  (with  amplitude  held  constant)  the  disturbance 
minimization  controller  of  this  section  did  reduce  the  amplitude  of  the  sinu¬ 
soidal  contribution  to  the  total  error,  as  compared  to  both  the  U4  ■  0  case 
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Figure  26.  Error  versus  time,  XSp-(10,0)T,  w=2SIN(3t). 
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Figure  34.  Percent  reduction  in  amplitude  of  sinusoidal  error  contribution 
alT  function  of  external  disturbance  amplitude,  a«=10  RAD/SEC, 
Xsn-(10,0)T. 
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Figure  36 .  Percent  reduction  in  amplitude  of  sinusoidal  error  contribution 
as  function  of  external  disturbance  amplitude,  a-3  RAD/SEC, 
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and  the  case  with  u<j  designed  to  minimize  ea8.  For  frequencies  below  about 
1.5  rad/sec,  the  controller  designed  to  minimize  e88  gives  better  performance. 
At  frequencies  above  about  25  rad/sec.  It  was  difficult  to  measure  performance 
accurately  due  to  the  small  amplitude  and  Increasingly  erratic  nature  of  the 
sinusoidal  component  of  the  response  produced  by  the  u<j  of  this  section. 

VIII.  DISTURBANCE  MINIMIZATION  WITH  w(t)  -  ce<*t 

This  section  will  apply  the  lsobasis  design  technique  to  the  state  set- 
point  regulator  example  with  an  external  disturbance  of  the  form 

w(t)  -  ce«t.  (84) 

For  this  case,  th»  a  priori  assumption  is  made  that  u^u  is  of  the  form 

u<jw(t)  “  bce0ct  -  bw(t) ,  (85) 

and  u^g  is  as  shown  in  Equation  (42). 

That  part  of  the  solution  of  Equation  (17)  for  e(t)  which  is  due  to  the 
last  term  on  the  right-hand  side  of  Equation  (17)  will  again  be  found,  with 
u<jw  as  given  by  Equation  (85).  In  this  case  then,  one  has 


/t  ~ 

eA(t-T)(F+Bb)ceai<iT 


-  -  (<zI-A)“lreaIt-eAt](F+Bb)c.  (86) 

If  only  the  response  after  transients  have  settled  out  is  considered.  Equation 
(86)  can  be  rewritten  as 


Ae^  -  -£“1(F+Bb)w  , 


where  A  ■  at  -  A. 

From  Equation  (87)  one  finds  that  that  b  which  is  itself  of  minimum 
norm  and  which  minimizes  the  norm  of  Ae„  is  given  by 

b*  -  -$-lBw)  tl-1-Fw.  (88) 


Evaluating  the  terms  in  Equation  (88)  results  in 


at  -  A 


’•“l. 

>.28j 


0.64T 


ll] 


ra-0.28 


a2  +  1.72a  +  3.28 


3a2-8.64crt-8.36 
5  a2-6  crt-5 


The  controller  defined  by  Equation  (85),  using  Equation  (93),  and  by 
Equation  (42)  was  added  to  the  set-point  regulator  digital  simulation. 

Results  are  presented  In  Figures  38  to  45  for  a  case  with  the  origin  as  the 
target  setpoint.  Data  Is  Included  on  each  figure  for  the  uj  ■  0  case  and  the 
case  with  u<j  to  minimize  e^.  From  Figures  38  to  40,  for  cases  with  negative 
exponents  In  w(t),  it  can  be  seen  that  the  lsobasls  controller  design  provides 
the  lowest  error.  From  Figure  41,  however,  It  can  be  seen  that  as  a  becomes 
Increasingly  negative,  the  lsobasls  controller  loses  Its  advantage.  Figures 
42  through  45  present  results  for  cases  with  positive  exponent  values  In  w. 

The  lsobasls  controller  provided  the  lowest  overall  error  In  all  cases.  For 
the  case  In  Figure  45,  the  controller  designed  to  minimize  eB8  produced  an 
error  twice  as  large  as  for  the  results  shown.  Figures  46  through  49  present 
results  for  two  non-zero  target  set-points.  The  results  are  mixed  for  these 
four  cases,  but  the  lsobasls  controller  can  provide  smaller  errors  In  some 
cases. 
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Figure  38.  Error  versus  time,  Xsp  *  ( 0,0 )T,  w(t)  -  lOe 
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IX.  SUMMARY  AND  CONCLUSIONS 


The  purpose  of  this  report  was  to  present  results  of  the  application  of  a 
disturbance  minimization  control  design  technique,  called  the  Isobasis  design 
technique,  to  a  state  set-point  regulator  problem  with  time-varying  external 
disturbance  inputs.  The  isobasis  design  technique  makes  the  a  priori  assump¬ 
tion  that  the  disturbance  control  vector  is  composed  of  some  combination  of 
the  same  basis  functions  which  describe  the  disturbances. 

Section  VI  presented  results  for  disturbance  minimization  when  an  exter¬ 
nal  disturbance  of  the  general  form  w(t)  -  c0+c^t  was  applied  to  the  plant 
described  in  Section  IV.  In  Part  A,  an  allocated  disturbance  minimizing 
control  vector  was  developed.  The  norm  minimization  and  critical-state 
variable  methods  were  illustrated.  It  was  shown  that  the  isobasis  technique 
provided  no  advantage  over  a  controller  designed  under  the  assumption  that  w 
was  a  constant.  In  Part  B,  an  unallocated  disturbance  minimizing  control  vec¬ 
tor  was  developed.  It  was  shown  that  with  this  approach,  the  isobasis 
designed  controller  can  provide  performance  equal  to,  or  better  than,  that 
provided  by  the  controller  designed  with  w  assumed  constant.  The  performance 
Improvement  does  not,  however,  occur  in  all  cases. 

Section  VII  presented  results  for  the  isobasis  technique  when  an  external 
disturbance  of  the  form  w(t)  ■  was  applied  to  the  plant.  It  was  shown 

that  the  lsobasis  designed  controller  can  provide  significant  reduction  in  the 
amplitude  of  the  sinusoidal  component  of  the  error.  Results  were  presented 
for  two  different  target  state  set-points  to  show  the  performance  of  the 
controller  as  a  function  of  the  frequency  of  w(t)  and  of  the  input  amplitude 
of  w(t). 

Section  VIII  presented  results  for  the  isobasis  technique  when  an  exter¬ 
nal  disturbance  of  the  form  w(t)  ■  ce  was  applied  to  the  plant.  It  was 
shown  that  the  isobasis  designed  controller  can  result  in  Improved  performance 
for  this  external  disturbance  also,  but  again,  not  consistently* 

As  these  example  results  have  shown,  the  isobasis  design  technique 
obviously  did  not  produce  the  best  disturbance  minimizing  controller  for  all 
external  disturbances  examined  or  even  for  all  target  set-points  with  the  same 
external  disturbance.  It  was  shown  that  it  does  have  the  potential  of  pro¬ 
ducing  a  controller  which  will  perform  well  in  reducing  the  error  between  the 
plant  state  and  the  target  state  set-point.  Its  application  should,  there¬ 
fore,  be  considered  in  cases  involving  time-varying  external  disturbances. 

In  the  Appendix,  a  definition  is  given  for  the  "utility"  U  of  time-varying 
external  disturbances.  Examples  are  presented  to  illustrate  the  concept  of 
disturbance  utility.  Two  external  disturbances  were  considered:  one  was  w(t)  * 
co  +  clt»  t*ie  ot^er  was  w(t)  “  csin(ott). 
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APPENDIX  A 


UTILITY  OF  TIME-VARYING  EXTERNAL  DISTURBANCES 


In  the  case  of  state  set-point  regulation/stabilization  problems 
involving  a  linear,  time-invariant  plant  and  a  constant  external  disturbance, 
it  was  shown  in  Reference  13  that  the  steady-state  error  is  the  residual  part 
of  the  sum  of  the  set-point  and  external  disturbance  vectors  which  lie  in  the 
orthogonal  complement  of  R(B)  and  is  hence  unabsorbable  by  the  control  vector. 
It  was  also  shown  that  the  steady-state  error  may  be  reduced,  independently  of 
the  control  vectors,  by  a  fortuitous  combination  of  the  set-point  and  external 
disturbance  vectors,  except  in  cases  where  the  set-point  is  the  origin.  The 
"utility"  U  of  the  external  disturbance  was  defined  in  Reference  13  as 

**"11  ^s  |  |  w«o  “  |  |  ess  |  |  wi*o*  (A-l) 

If  U  >  0,  then  w  can  actually  aid  in  reducing  esa.  On  the  other  hand,  if  U  < 
0,  then  the  presence  of  w  increases  the  value  of  e8a.  Some  examples  of 
disturbance  utility  were  presented  in  References  13  and  14. 

In  the  case  of  time-varying  external  disturbances,  as  considered  in  this 
report,  the  definition  of  U  given  by  Equation  (A-l)  must  be  modified  to 

0  -  II  -  ||e<t>  ||U)lo.  (A-2) 

Depending  upon  the  system  and  the  particular  w(t),  it  is  possible  that  U  could 
be  always  positive,  always  negative,  or  positive  until  a  certain  time  and  then 
negative  thereafter,  or  vice  versa.  Also,  the  utility  could  exhibit  a 
periodic  behavior  as  a  function  of  time,  being  alternately  positive  and 
negative.  Several  examples  will  be  presented  in  this  section  to  illustrate 
disturbance  utility  for  the  plant  model  given  by  Equation  (18)  and  several 
types  of  time-varying  external  disturbance. 


For  the  first  example,  the  external  disturbance  will  be  given  by 
w(t)  -  1  +  t.  (A-3) 

From  Equation  (17),  with  u<j  ■  0,  the  solution  for  e(t)  is  found  as 


e(t) 


eAte(0)  - 


/t  ~ 

eA(t-x)[Axgp  +  Fw(-c)]dx 


(A-4) 


Let  the  target  set-point  be  xsp  ■  (-5,  0)T.  With  w(t)  *  0,  one  has 

II  *t)  llw-o  "  ||  /^V^^H-H-A-lAXgpll  (A-5) 

o 

if  the  initial  transients  given  by  ;*Ate(0)  are  ignored.  Making  the 
appropriate  substitutions  into  Equation  (A-5),  the  result  is 


|  £(t)  |  |  w.0  -  | | (0.427,  9.145)?  |  |  -  9.154. 


(A-6) 


With  w  as  given  by  (A-3),  the  norm  of  e(t)  is  found  as 


(A-7) 


|  |  ***)  |  |  wrfo  "  |  |  r  eXCt-cJIAxap  +  F(l+T)]d<  |  | 

|  j -A-1(Axsp  +  F)  -  A”1 (A-1  +  It)F  |  | 

-  V 1.5t2  -  21. 5t  +  79.4  . 

In  order  to  find  the  time  at  which  U  changes  sign  one  sets 

1 1  I  I  w-o  ”  1 1  E(f)  1 1  w^o  "  0.  (A-8) 

Making  the  substitutions  from  Equations  (A-6)  and  A-7)  into  (A-8)  results  in 

1.5t2  -  21. 5t  -  4.4  -  0,  (A-9) 

which,  when  solved,  yields  a  positive  time  of  t  ■  14.5  seconds.  In  this  case 
then,  at  t  ■  14.5  seconds  the  utility  of  the  external  disturbance  goes  from 
positive  to  negative,  l.e.,  from  aiding  in  error  reduction  to  increasing  the 
error  magnitude.  A  digital  simulation  was  written  for  this  example  and  the 
result  is  shown  on  Figure  A-l.  As  can  be  seen,  e(t)  for  the  case  with  w  -  1+t 
exceeds  e(t)  for  the  case  with  w  ■  0  after  t  *  14.5  seconds. 

If  the  target  set-point  is  assumed  to  be  xSp  -  (10,  0)^  then 

)  |  eCt)  |  |  w-o  “  |  | (-0.854,  -18.29)T  ||  -  183.1  (A-10) 

|  |-«Ct)  Hw*,  -  1 1  (0.1097t  -  0.497,  -1.2 19t  -  18.566)T  || 

-  Vl.5t2  +  45.16t  +  345.1  .  (A-ll) 

From  Equation  (A-8)  one  finds  that  for  this  set-point, 

1.5t2  +  45.16t  +  9.9  -  0,  (A-12) 

and  there  are  no  positive  values  of  t  which  will  satisfy  Equation  (A-12).  For 
this  set-point  then,  the  external  disturbance  always  acts  to  increase  the 
error  (since  at  t  ■  0,  0  ■  -9.9).  This  case  was  also  simulated  and  the  result 
is  shown  on  Figure  A-2.  As  can  be  seen,  e(t)  for  the  case  with  w(t)  ■  1  +  t 
is  never  smaller  than  e(t)  for  the  case  with  w  -  0. 

As  another  example,  let 

w(t)  ■  sin(t)  (A-13) 

with  x8p  -  (10,  0)T.  In  this  case,  one  would  have 

||  e(t)  Uw^o  -  ||“  f  e^(t“T)d'c(AxSp)  -  J'  eA(t“'t^simdT(F)  || 

o  o 

■  |  |  A'^Axsp  -  ^-1[-Asin(t)  -  Icos(t)]F  | | , 


(A-14) 
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Figure  A-l. 


Error  versus  time  for  cases  with  w»0  and  w»l+t 
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where  A  ■  +  c?I.  When  Equation  (A-14)  is  evaluated,  the  result  is 

II  /0.854  -  0.355sin(t)  -  0.232cos(t)\  || 


e( c  )  w^o 


18.79  +  l.llsin(t)  -  1.336cos(t) 


(A-15) 


The  norm  of  the  error  when  w  -  0  is  given  by  Equation  (A-10).  The  utility  of 
the  external  disturbance  of  Equation  (A-13)  and  the  given  set-point  is  found 
from  Equations  (A-15)  and  (A-10)  to  be 


U  -  183.1  - 


0.854  -  0.355sin(t)  -  0.232cos(t) 
18.79  +  l.llsin(t)  -  1.336cos(t) 


(A-16) 


The  digital  simulation  used  to  generate  the  data  in  Figures  A-l  and  A-2 
was  modified  to  include  the  w(t)  as  given  by  Equation  (A-13)  and  to  calculate 
U  as  given  by  Equation  (A-16).  Figure  A-3  presents  U  versus  time  and  Figure 
A-4  presents  the  corresponding  error  versus  time  data  for  this  second  example. 
If  Figures  A-3  and  A-4  are  overlaid,  the  agreement  can  be  seen  between  the 
regions  of  positive  utility  and  the  regions  where  the  error  with  w  j4  0  is  less 
than  that  with  w  ■  0.  Figures  A-5  and  A-6  present  similar  data  for  a  case 
with  w  ■  sin(t)  and  xgp  ■  (-5,  0)T.  Again,  it  can  be  seen  that  the  utility  of 
w(t)  alternates  between  positive  and  negative. 

Figure  A-7  is  a  plot  of  the  state  space  showing  the  external  and  set- 
point  disturbance  vectors  and  the  line  of  action,  R(B),  of  the  control 
corresponding  to  the  case  of  Figure  A-5.  Note  that  the  vector  representing  Fw 
has  an  arrowhead  at  each  end.  This  represents  the  fact  that  w(t)  *  sin(t)  and 
w(t)  varies  between  +1  and  -1  as  a  function  c£  time.  The  component  of 
Axap  lying  in  R(B)-1-  is  shown  and  is  denoted  by  a.  Since  w(t)  varies  with 
time,  the  component  of  Fw  lying  in  R(B)X  will  also  vary.  This  component  is 
also  shown  in  Figure  A-6  and  is  denoted  by  f. 

In  Reference  13,  two  criteria  were  given  which  must  be  satisfied  in  order 
for  w(t)  to  exhibit  a  positive  utility.  The  first  criterion  Is  a  magnitude 
criterion. 


<  2TT  , 


and  the  second  is  an  angle  criterion, 
90°  <  0  <  270°  , 


(A-l 7) 


(A-l 8) 


where  0  is  the  angle  between  the  two  vectors  (2a+f)  and  (f).  From  Figure  A-7, 
it  can  be  seen  that  Equation  (A-17)  is  satisfied  for  all  t.  The  angle 
criterion,  however,  is  only  satisfied  when  vw  is  positive,  i.e.,  when  sin(t) 
is  positive.  Since  a  is  given  by  [13,14] 


*  ■  (I  -  BB^AXgp  -  (-4,  2)T  , 
the  magnitude  of  "S  is  found  to  be 

INI-4-47 


(A-19) 
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Figure  A-4.  Error  versus  time  for  cases  with  wO,  w=SIN(t) 
with  XSp»(10.  f  0. )  iJ. 
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If: 


From  the  magnitude  criterion  Equation  (A-17),  even  when  9  satisfies  the  angle 
criterion,  If  ||f  ||is  greater  than  |  |  2a  ||,  the  utility  of  w(t)  will  be 
negative* 

In  order  to  Illustrate  how  utility  changes  with  w(t),  two  additional  runs 
were  made  for  the  xg^  ■  (-5,  0)T  case*  Figure  A-8  shows  the  results  for  w(t) 

■  10sin(t),  i.e.,  ||f  ||*  4.47,  and  Figure  A-9  shows  the  results  fo£  w(t)  ■ 
20sin(t),  i.e.,  ||7  ||“  8.94.  As  can  be  seen,  as  the  magnitude  of  f 
approaches  the  limit  set  by  Equation  (A-17)  the  utility  of  w(t)  approaches  a 
condition  where  it  will  be  zero  for  all  t  (since  the  error  with  w  +  0  is 
approaching  a  condition  where  it  will  always  be  greater  than  for  the  case  with 
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